Pairing, Ferromagnetism, and Condensation of a normal spin-1 Bose gas 
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We theoretically study the stability of a normal, spin disordered, homogenous spin-1 Bose gas 
against ferromagnetism, pairing, and condensation through a Random Phase Approximation which 
includes exchange (RPA-X). Repulsive spin- independent interactions stabilize the normal state 
against both ferromagnetism and pairing, and for typical interaction strengths leads to a direct 
transition from an unordered normal state to a fully ordered single particle condensate. Atoms with 
much larger spin-dependent interaction may experience a transition to a ferromagnetic normal state 
or a paired superfiuid, but, within the RPA-X, there is no instability towards a normal state with 
spontaneous nematic order. We analyze the role of the quadratic Zeeman effect and finite system 
size. 



Introduction. — The interplay of superconductiv- 
ity/superfluidity and magnetism is fundamental [H. Ex- 
periments in ultra-cold spin-1 gases Q-lll have begun to 
explore this physics, elucidating the subtle connections 
between Bose condensation of single particles and com- 
peting/complementary orders such as pair condensation, 
ferromagnetism, and liquid-crystal like nematicity 
This has become a model system for thinking about ex- 
otic spin textures and topological defects |9|; 12i. and the 
dynamics of quantum phase transitions |13l . Il4 |. How- 
ever, the finite temperature 3D phase diagram is a mys- 
tery, with earlier works producing contradictory results 
15l4l7l|. Here we clarify the situation by using a well- 
controlled approximation [the random phase approxima- 
tion with exchange (RPA-X)] to calculate the instabilities 
of the normal state. 

A particularly dramatic feature of the spin-1 Bose gas 
is that it sup ports a bosonic analog of the BCS tran- 
sition |18l42l| . Somewhat counterintuitively, the paired 
state is less ordered than a single particle condensate, and 
is found when both the spin independent and dependent 
interactions are repulsive. In addition to its theoretical 
importance, this feature makes it an interesting paradigm 
to keep in mind while exploring the mechanisms for su- 
perconductivity in systems such as high-Tc cuprates, C- 
60 and polyacenes where the interactions are believed to 
be repulsive [j^ . 

The Hamiltonian of a spin-1 Bose gas is the sum of a 
kinetic and interaction term, H = Hkin +'Hint- In the 
presence of a magnetic field in the z direction, the kinetic 
term has the form = J^ka ^kaa\„aka, where is 
the annihilation operator for a boson with momentum 
k and spin projection a ~ —1,0,1. The dispersion is 
ffeo = k'^/2m — /i, ek±i = /2m — fi + q zL p, where p/q 
are linear/quadratic in the magnetic field. There is no 
spin-orbit coupling, allowing us to eliminate the linear 
Zeeman effect (p) by working in a rotating frame. Off- 
resonant microwave light allows the quadratic Zeeman 
field q to be tuned, taking on positive and negative values 
[2^. Assuming short range interactions, symmetry forces 



TABLE I: Orders in spin-1 gas. 
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where the greek indices denote the spin projection and 
i^a{r) = y J2k ^^''^^ka is the the boson field operator. 

The two coupling constants, cq and C2 represent spin 
independent and spin dependent interactions The S op- 
erators denote 3x3 spin-1 matrices. The interac- 
tions are expressed in terms of the microscopic scat- 
tering lengths in the spin-0 (oq) and spin-2 (02) chan- 
nels and atomic mass m as: cq = 47r(ao + 2a2)/3m and 
C2 = 47r(a2 — ao)/3m. 

Two atoms are typically used in these experiments: 
^''Rb (c2 < 0) and ^■^Na (c2 > 0). In all experiments so 
far, spin independent interactions are repulsive (cq > 0) 
and Co S> \c2\. We find that the phase diagram is fea- 
tureless in this regime, motivating us to study the more 
general case where the interactions are comparable in 
magnitude. Perhaps ^Li, rare earth atoms, or alkali-earth 
atoms will have scattering parameters in this regime. Al- 
though dipolar interactions are believed to play an impor- 
tant role in the low temperature quasi-2D experiments of 
Vengalattore et al. Q , we neglect them here, as they are 
much too weak to influence the stability of the normal 
state. 

The spin-1 gas can present several types of order, 
summarized in Table HI Single particle condensate or- 
der (C) is always accompanied by either ferromag- 
netic (FC) or nematic order (NC). Mean field exam- 
ples of these condensed states are \FC) = (?/'|)^|0), 
and \NC) = (V'J)^|0), with the former seen in ^'^Rb, 
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FIG. 1: 3D Phase diagram at q = within the RPA-X: 

Thick solid/dashed lines are the ferromagnetic and pairing 
transition temperatures measured from the ideal Bose gas 
transition ksTBEC = If ( 3c(3/2) 
as a function of spin-dependent interaction C2 



'^''^ and scaled by n^''^ao 



Thin solid 

line shows the instability towards single particle condensa- 
tion (ferromagnetic or polar depending on the sign of C2). 
For C2 > 0.5co, the normal state is unstable to a rotation- 
ally symmetric paired singlet phase (P) with a Tc > Tbec- 
For C2 < —Co/3, the normal state becomes unstable to a fer- 
romagnet (F). For |c2| smaller than these threshold values, 
there is a direct transition from a disordered normal state 
to a ferromagnetic condensate (FC) or a nematic condensate 
(NC) . Tick marks on the upper frame illustrate the scattering 
lengths for ^^'^Rb and ^^Na. 



and the latter in ^'^Na. The singlet state from [2J| 
with all particles in fc = 0, \S) = ((V'o k=o)^ ~ 
^''^1 k=o^-i k=o)^^^l'^)' off-diagonal single-particle 

order [ie. lini|r_r'|^oo('*/'o(^)V'o(''')) 7^ 0] in the 

thermodynamic limit should be considered as an NC 
state [25|. An example of a P state would be a con- 
densate of small singlet pairs \P) = K^^^\0}, where 

K = J2k ('^ok'^o-k ^ ■ Unlike \S), the state 

|P) has no off-diagonal single particle order. Paired 



states for which 



^Ok"o-k 



2a 



lk"-l-k/ 



7^ posses ne- 



matic order. We find no instabilities towards paired 
states with ferromagnetic order. 

Although the 2D phase diagram is well established 
(with an algebraically ordered P state at any finite tem- 
perature when q = 0) [lol . 16 1, contradictory results have 



appeared concerning the 3D phase diagram. Both Gu 
and Klemm [l^ and Kun Yang [3| erroneously found 
that arbitrarily weak attractive interactions drive a fer- 
romagnetic instability with Tf > Tbec- Kis-Szabo, 
Szefalusy and Szirmai gave a more thorough argu- 
ment, finding a finite threshold for this instability. We 
extend their calculation to incorporate exchange physics, 
including C2 > 0, C2 < 0, and the quadratic Zeeman ef- 
fect. 

Formalism. — We calculate the longitudinal (xz) and 
transverse spin (x±) and pairing (11) susceptibilities of 
the homogeneous interacting spinor Bose gas using a 
Hartree-Fock Random Phase Approximation (RPA-X). 
A divergence of the zero frequency, long wavelength sus- 
ceptibility, (k = 0, w = 0) =0 signals an instability in 
that channel. Within the RPA-X there is never an insta- 



bility in the nematic channel which is not simultaneously 
accompanied by single particle or pair condensation. 
The relevant response functions are 

Xl'p(*) = 7(E akW«7k+p(i)«k(0)«"q-p(0)) (2) 

k.q 

= 7(E4kW4p_kW«^^q(0)«cp-q(0)) (3) 

where t > 0, and the greek subscripts denote spin in- 
dices and p is the momentum 2^ 27 1. The longitudinal 
and transverse spin correlation functions are Xzp(t) = 
-i(^.p(t)^._p(0)) - -Ea.5=±i(-l)"'xiipW, and 
X± = -z(5+p(t)^-p(0)) = xj? + X^o-l + X"o-; + X'll 
where Sf,p{t) = Eq «p+q/2(0'^^ap_q/2(i), and = 
{z,±}. 

In the RPA-X, the susceptibility of the interacting gas 
is determined from the non-interacting susceptibility by 
summing over all repeated direct, and exchange interac- 
tions. 

ix'"'^^)!} = ix')Zs^,Sp, + E(/);^?v;^i!(x«^^)r^(4) 
(n«pA).. ^ (n«)^^,<5„,<5,, + 5](n«);5!;v;i'j(n«^^)i;:;(5) 



The interaction potential V^J' of Eq. ([T]) , which includes 
both direct and exchange graphs, is explicitly given in 
the supplementary material. 

The non-interacting Green's functions are diagonal in 



spin space: (x )l1p{t) 
•q ~ a and 7 = /3, 

ix'&P,^) = 



0, and {W)ll^{t) 



unless 



(27r)3 uj ~ {ek+pp - eka) 
d^k 7i(£k,«) + »-(ek+p,/3) 
(27r)3 CO - (ek+p/3 + Cka) 



(6) 
(7) 



Here n{eka) = [e^'^'"' — 1)"^ is the Bose-Einstein distribu- 
tion at temperature T = For a non-interacting gas, 
the spin susceptibility x°, pairing susceptibility 11° and 
compressibility all diverge as /i — > from below, marking 
Bose-Einstein condensation. 

At fc, a; = 0, these non-interacting response func- 
tions may be written in terms of the polyloga- 



rithm functions gy{z) 



0,w 



4irA: 



0) 



^gi/2(e^(^-'?)), (x°)J?(0,0) 



:(T/g)[.93/2(e''(^-')) - 33/2(6'^^)], n^;(0,0) = 



H — q for a ~ ±1 and 
= /3 = 0. The thermal 



-^5i/2(e^^"«), where ^eff = 
j3 — =pl, and /ioff = for a 
wavelength is Ay = ^J2'K /rakBT . The calculations are 
detailed in the supplementary material. 

To detect fcrromagnctism we consider the response 
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functions (sec supplementary materials) 

2(x°)_r}(k,c.) 



(a) 



1.13 



xr^(k,o) 
xr^(k,o) 



l-(co + 3c2)(x°)_rj(k,0) 

2(x°)}g(k,0) 
l-(co + 3c2)(x")J?(k,0) 



(8) 
(9) 



To detect pairing, it suffices to consider the singlet pair- 
ing susceptibility, 8 ^ (nSj[j - 2n_i"J)^^^, 



e = 



n+ - 2n" + n"n+(co - 202) 



1 - (co - C2)n+ - cqHO + (Co - c2)(co + c2)n+no ■ 

with = [U")^^ and n+ = (nO)_J~J . when g = 
(nO)gg = {U")_l^l = (n)0 and this expression simplifies 
to cx f - (co - 2c2)(n)0. 

Results. — Repulsive spin independent interactions (cg) 
suppress both ferromagnetism and pairing in Bose sys- 
tems. This should be contrasted to fermions, where due 
to the opposite sign of the exchange term, repulsive inter- 
actions enhance ferromagnetism, giving rise to the Stoner 
instability [28| , even in the absence of any spin dependent 
interactions. 

From Eq. HI and the fact that x°(0, 0) < we see that 
the spin susceptibility only diverges when C2 < with 
|c2| > Co/3. Similarly, at </ = 0, the pairing susceptibility 
only diverges when C2 > with C2 > co/2. For weak 
interactions (|c2n| ksT), these instabilities occur near 
^ = 0. Expanding the susceptibilities for small /i at g = 
gives that, to leading order, the magnetic instability 
occurs at 

W = ^^r^ - 4.84 (i - nV3„„ (10) 

Vir = -^-^ =6.44 --- n^/^ao. 

Jbec \Co ^/ 

Hypothetically, taking n = lO^^cm-"^, uq = IOOgb, and 
|c2| ^ Co, we find Tc — Tbec ^ lOnK. The g = phase 
diagram is summarized in Fig. [TJ 

We now explore the role of the quadratic Zeeman ef- 
fect: g < favors magnetism in the ±z direction (Fy 
- Ising order - signalled by diverging Xz), and pairs in 

the mp = ±1 states (NP_l : 2|(VIV'li)l > l(V'JV'J)l); 
while g > favors magnetism in the x — y plane 
(Fj^ - x-y order ~ diverging x±): and mp = O pairs 
(NP|| : KV-JV^i)! > 2|(V'IV'li)l)- Finite q also shifts 
the BEC transition temperature: the density is given by 
nk^ = 33/2(6^'') -I- 2p3/2(e^'^''-«)), with condensation at 
/i = g for g < and at = for q> {). For small q one 
finds 



T^tl=ntc+^^ntc^ (11) 

with C = 0.3 for g < and ^ = 0.6 for g > 0. 
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FIG. 2; Instability with C2 < 0: (a): Instabilities of the 
unordered normal state with C2 = — cq as a function of q. 
Solid curves give the for a non-condensed ferromagnetic 
gas, normalized to TBEc\q=o (defined in Fig[l] caption). At 
some lower temperature, one expects a transition to a ferro- 
magnetic condensate (FC). For q < 0, this Tc always exceeds 
the ideal Bose gas transition temperature. For q> Q, the ideal 
gas temperature meets the Tc for ferromagnetism at some fi- 
nite q (marked by x). Beyond this point, the normal state is 
unstable to forming a polar condensate, (b): Location of x 
as a function of interaction strength. 



Figure [2l^a) illustrates the phase diagram for C2 < 
0; where the only relevant instabilities are ferromag- 
netism and single particle condensation. For g < 
and |c2| > Co/3 an Ising ferromagnetic instability always 
precedes condensation. For g > there is a thresh- 
old g below which x — y ferromagnetism precedes con- 
densation. The dependance of this threshold on C2 is 
shown in Figure [Hb). For C2 near — co/3, one finds: 
g^ = r«=° (10.6(aon)i/3a)2, where a = 1/3 - |a2|/ao. 

Figure [3Ja) illustrates the phase diagram for C2 > 0, 
where the only relevant instabilities are pairing and single 
particle condensation. Finite g enhances single particle 
condensation, and for a given g, there is a threshold value 
of C2 required to find a pairing instability. For g > this 
threshold becomes arbitrarily large as g — > 00, but for 
g < one always has a pairing transition if C2 > cq. 
Setting ^ = 0, and taking the limit (n°)_J~J — )■ 00 for 



g < and (n°)oo — > 00 for g > 0, we calculate these 
threshold values (Fig. [3^b)): 



ge = 12.74T«^"c("'^'ao)'(l-22;)(l + .x) 



1, '7>0 I 
9<0/ 
(12) 

where x = C2 /co . 

Experimentally, the states discussed in FigjS] may be 
somewhat distinguished by the fact that both the con- 
densed phase NCx for large g < 0, and the paired phase 
NP_L, have Us = iii + n_i — 2no > 0. In the singlet pair, 
this quantity is identically zero. Studying momentum 
distributions can distinguish between the single particle 
and paired condensates. 

Finite Size effects. — Similar to p]| . we can estimate 
the role of finite size effects by looking at the instabilities 
at finite k = 27r/L, where L is the size of the cloud. 
These finite size effects arc crucial in the scalar gas with 
attractive interactions, where there is no Bose-Einstein 
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FIG. 3: Instabilities with C2 > 0: (a): Instabilities of the 
unordered state with C2 = 0.75co. At large negative q, the 
normal state is unstable towards a single-particle condensate 
(NCx) with a spinor order parameter ^ — {e"^,0, 1}), with 
arbitrary 0. For — l^ciKx) < q < 15021(0) the instability is 
towards pairing. At q — the paired phase is a spin singlet 
with no spin fluctuations. For g 7^ 0, the paired phase has spin 
fluctuations in the x — y plane. For large q > 0, the normal 
state is unstable towards a polar condensate with spinor ^ — 
{0, 1,0}. Dotted lines are the ideal gas condensation Tc as a 
function of q. (b): Details of dashed rectangle in (a), showing 
instabilities towards pairing, (c): Critical values of g at which 
pair instability gives way to single particle instability (i.e. 
locations of x and o) for different interaction strengths. 



condensation in the thermodynamic limit. Here we find 
only small corrections to location of the threshold value 
of |c2|/co. To leading order in kAr <C 1, the threshold 



for ferromagnetism is 02/00 = 



kA 



1 _ 

'3 



J^BEQ Tff\y[\Q the 

127r^ao ' 



threshold for pairing is 02/00 — 2 < sn-^ao 
are no greater than 10% for a system of size L ~ lOOyum. 

Conclusions. — The presence of competing magnetic 
and off-diagonal long range orders in a spin 1 gas pro- 
duces an extremely interesting phase diagram with fer- 
romagnetic, nematic and paired phases. Using the RPA- 
X, we have quantitatively studied the instabilities of the 
normal state, identifying the temperatures and interac- 
tion strengths at which the disordered normal state be- 
comes unstable to a symmetry broken phase. We find 
that the finite temperature phase diagram is featureless 
unless the interaction strengths governing spin (02) and 
charge physics (co) are comparable in magnitude. This 
is due to the exchange enhancement of identical particle 
scattering. 

A number of probes can be used to distinguish the 
nematic phases and detect pairing. These include optical 
birefringence (29j . momentum distributions via time-of- 
flight, noise correlations [1^, and the nature of vortices. 

Finally, we remark that the key bottleneck to realizing 
this interesting phase diagram is that of finding atoms 
with spin independent and spin-dependent interactions 
of comparable magnitude. The search for such atoms is 
an active area of research, and we hope that our work 
motivates this effort. 
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SUPPLEMENTARY MATERIALS FOR 
"INSTABILITY OF A SPIN-1 BOSE GAS TO 
FERROMAGNETISM AND PAIRING" 

Analytic structure of non-interacting response 
functions. — Wc develop the analytie structure of Xq^" 
defined in Eql6]of the main text as: 



Rewriting the logarithm as an integral we get: 



dxli (k) 



(P'k n{ek,a) - n{ek+p,i3) 



(A-13) 



(ix/o(fc) 
(A-16) 

where /i/o(fc) ^ C^dk^ where Ai = 

—/3{fJ. — q) and Ao — -~/3fJ,. The analytic structure of the 
integral / has been extensively developed by Szepfalusy 
and Kondor who show that the integral can be writ- 
ten as an asymptotic series for long wave-lengths. Re- 
taining only the lowest order terms we find that the static 
response yields: 



For q ~ this simplifies to a constant xip^^) = 
f "('^'°)~"('^''+p) whose structure has been thor- 
oughly explored in [l| . The result to linear order in pAx 



m kgT 



(ff3/2(e^(^-''')-53/2(e^^)) 

(A-17) 



The non-interacting pair response is defined as: 



ITT ( 

log 



pA 



where gi,{z) is the polylogarithm function. For q ^ 0, 
first note that xf;(0,0) = x(0,0)U^p-9 for = 
{1,-1}. The transverse spin response however given by 
(x*')oi(Pj '-^) requires some work. One proceeds as follows: 
Integrating out the angular variables one finds 



Once again, the pair susceptibility for the scalar gas 
has been considered in 2l|. The result to linear order in 
/sAt is : 



n"(p,o) 



2i7r , 
+ — log 



-:^{9i/2{e^n^VWm\ (A-19) 
'(l + ^)f^-^/^\ 



pAr 



(1 



(X°)^?p(-) = 



ttA^P 7o 



dfcn(ej,J log 



+ '^(efco) lo 



2_ + fc 



1 Mr 



z_ + k 

Z-\- — k 
.z+ + k 



r (A-15) 



where z± — , 

2y/epkBT 

k as an expansion parameter 



± and fc fcAT/40F. We use 



It is easy to show that (n°)_i[~i[(0, 0) = 

nO(o,o)U^^_, and (n")g[] = n". 

Spin response in the RPA . — The spin response in 
the RPA is given by solving for the polarization tensor 
defined in Eq. [2] using Eq. [5] The interaction matrix V 
encompasses all direct and exchange diagrams and takes 
the form: 



/ 2(co + C2) 











Co + C2 











Co 


-C2\ 











Co + C2 











2C2 

























Co - C2 













Co + 


C2 








2C2 













Co + C2 











2co 











Co 


+ C2 











2C2 








Co 


+ C2 













Co - 


-C2 






















2C2 











Co + C2 













V Co - C2 











Co + C2 








2(co^ 


-C2)/ 



r 



From the full polarization tensor, one extracts the lon- 



gituginal and transverse spin susceptibility on which our 
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calculations are based. 

We now turn to the details of the pair response calcu- 
lation. 

Pairing response in the RPA. — The RPA response is 
given by Eq. [5] where ^ is a symmetric 9x9 matrix. 
However since pairing only occurs in the Sz = channel, 
it sufBces to consider the following subsystem 

n n^-i n^-i \ 
nt; no« n_?? (A-21) 

^^-11 00 1-1 / 

which is related to the non-interacting response 

u_l-l o\ 

nj][] via Eq.[5]and the 3 X 3 interac- 

u-uj 

(Co - C2 C2 \ 

C2 Co ■ 
C2 Co - C2 / 

Note that n_i-J = n-j_J and the physical response 
of the system to adding a ±1 pair is given by n_|[-J + 
n-j_J. The 3x3 system can be further reduced to a 
2x2 system of equations which then yields the singlet 
RPA response. 



Two limiting cases are worth considering. The first 
is at /i = q for q < when the non-interacting BEC 
transition occurs for the ±1 atoms. At these values the 
non-interacting response function n_j-J diverges, and: 

(X -Co + C2 + (co - 2C2)(C0 + C2)n0° (A-22) 

The second is at /i = which corresponds to the BEC 
transition for the atoms. At this point non-interacting 
response function IIoo diverges and: 

cx -co + (co - 2c2)(co + C2)n_J-i (A-23) 

Setting the R.H.S. of Eqs. (|X^r03)) to zero, using 
EqfTT] along with the functional forms of the response 
functions, yields the threshold value of C2 at which paired 
states result for any q. 



Electronic address: ls sn8@cornell.edul 
[1] E. J. Mueller and G. Baym, Phys. Rev. A 62 053605 
(2000). 

[2] P. Szepfalusy and I. Kondor, Ann. Phys. 82 1 (1974). 



